Chaplygin gas in light of recent integrated Sachs— Wolfe effect data 
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We investigate the possibility of constraining Chaplygin dark energy models with current Inte- 
grated Sachs Wolfe effect data. In the case of a flat universe we found that generalized Chaplygin 
gas models must have an energy density such that fl c > 0.55 and an equation of state w < —0.6 at 
95% c.L. We also investigate the recently proposed Silent Chaplygin models, constraining Q c > 0.55 
and w < —0.65 at 95% c.L. Better measurements of the CMB-LSS correlation will be possible with 
the next generation of deep redshift surveys. This will provide independent and complementary 
constraints on unified dark energy models such as the Chaplygin gas. 



PACS numbers: 

I. INTRODUCTION 

A decade of impressive progress in cosmology has fi- 
nally led to the so-called ACDM model of structure for- 
mation, based on inflation, cold dark matter and a cos- 
mological constant (see e.g. QJ). However, this can't be 
seen as the ultimate result, since its features include a 
massive amount of energy of unknown origin in the uni- 
verse, usually split in two categories called dark matter 
(DM) and dark energy (DE) (see e.g. 0). 

In fact, there are strong experimental constraints that 
have induced the introduction of new items in the cos- 
mic energy balance, and these constraints arise from two 
different classes of observations: the ones based on the 
study of the property of the matter and its clustering 
have suggested the presence of DM, while the study of the 
cosmic microwave background (CMB) and supernovae la 
(SNaela) has led to a wide consensus about DE. 

Nevertheless, our theoretical knowdledge of the dark 
universe is still quite poor since it is not directly detected 
but only inferred from its cosmological implications, and 
many models for DM and DE are not yet ruled out by the 
observations. For the sake of simplicity, one could then 
make the ansatz that these two entities are in reality 
just two distinct manifestations of one thing. This is the 
purpose for which have been introduced the unified dark 
energy models (UDE), also called quartessence. 

An alternative to quintessence model would be an ex- 
otic fluid capable to develop a negative pressure at late 
times, while approximating the dust matter behaviour 
earlier on, with a smooth transition: these are exactly 
the features of the Chaplygin gas, a fluid introduced in 
1904 Q for aerodynamics, that presents some interesting 
properties recently discovered for cosmology Q. 

The Chaplygin gas is characterized by an equation of 



*Electronic address: tommaso.giannantonio@port.ac.uk 
^Electronic address: alcssandro.mclchiorri@romal.infn.it 



state of the form 



A 

Px = — 
Px 



with a = 1, while a generalized Chaplygin gas has < 
a < 1 tii an d A is a constant with dimensions M^ 1+a \ 

The Jeans lenght of the GCG perturbations is first 
similar to the matter one, and then the instability is re- 
moved when the behaviour approximates a cosmological 
constant; this produces a large integrated Sachs- Wolfe 
(ISW) effect jfj on the CMB anisotropies, significatively 
different from the one expected in ACDM models. Re- 
cent combined analysis of CMB anisotropies and Large 
Scale structure data have severly constrained the Chap- 
lygin gas model (see for example IHj). However, several 
cosmological backgrounds may mimick an enhancement 
in the large angular scale anisotropy as the one produced 
by the ISW like, just to name a few, gravitational waves 
or cosmic strings. It is therefore important to verify this 
result by using independent and complementary datasets. 

A powerful and model-independent way to extract the 
ISW signal from large scale CMB anisotropy has been 
proposed in [9j by cross-correlating CMB maps with 
galaxy surveys. With the advent of the new WMAP re- 
sults, several cross-correlation analysis have been made 
with at least five detections at > 2ct level. It is therefore 
timely to investigate the impact of those measurements 
on the Chaplygin gas model. In the following of the pa- 
per, we use the recent detections of the ISW effect in the 
CMB anisotropies resumed in 0. 

In the next section we will discuss the Chaplygin gas 
model and the expected ISW signal. In section 3 we will 
analyse the current data and produce new independent 
constraints on the model and finally in section 4 we will 
derive our conclusions. 
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II. THE CHAPLYGIN GAS AND THE 
INTEGRATED SACHS-WOLFE EFFECT 

From equation ^ an d energy conservation it follows 
that for the Chaplygin (see e.g. [Tof') 



Px 



A 



D 



a 3(a+l) 



(2) 



where B is a constant with the same dimensions of A] 
this means that the equation of state for the GCG is of 
the form 



w(a) = 

that at present time is 



Px 
Px 



A 



pT' 



A 



A + B 



(3) 



(4) 



Using the physical parameters Qx and w instead of A and 
B, we can recover the expression of the evolving equation 
of state: 



w(a) 



1 



wa 3(a+l) 



- 1 



(5) 



From the equation above it is clear that at early times 
the fluid behaves as non relativistic matter while at late 
times it behaves as a fluid with equation of state w. 
The corresponding density evolution equation is 



p(a) 



1 



w 



a 3(a+a) 



(6) 



In the synchronous gauge the evolution equations for 
the density and velocity divergence perturbations, S and 
9, in Fourier space for the Chaplygin gas are (see e.g. 

D3) 



-(1 + w) [-0 + - ) + 3Hw(l + a)S 



6= -(l + 3wa)H9 



1 



-k 2 S 



(7) 



(8) 



where the derivatives are respect to the conformal time 
(d/dr) and aH — da/dr. As showed in below a char- 
acteristic scale k 2 > H 2 /(\aw\) S has oscillating (grow- 
ing) solutions for a > (a < 0). This behaviour makes 
the Chaplygin gas at odds with the current cosmologi- 
cal data when considered as unified dark matter model. 
The generalized Chaplygin gas can therefore be consid- 
ered only as a candidate for the Dark Energy compo- 
nent. Recently, Amendola et al. 01 have introduced a 



new version of the Chaplygin model, whith an additional 
entropy perturbation component T such that: 



Sp dp 



c 2 = 



(9) 



i.e. the effective sound speed of the cosmic fluid van- 
ishes, assuring clustering on small scales. In this case, 
the perturbation equations change to 



-(1 + w) + J) +mw8 (10) 



6 = -(l + 3wa)W (11) 

The Silent Chaplygin model is in good agreement with 
the current status of observations [TtJ ■ 

Differences between Chaplygin models may arise at 
late times, when their energy contribution becomes im- 
portant; for this reason, it is useful to study the con- 
sequences of both models on a typically late time phe- 
nomenon as the integrated Sachs- Wolfe effect 

When CMB photons pass through potential wells they 
can acquire a red- or blue-shift if the gravitational poten- 
tial is not constant in time; this is exactly what happens 
when in the universe energy density balance dark en- 
ergy becomes to be important. This effect can only arise 
at late times, when z < 4 (and fi m significantly differ- 
ent from one), and when galactic structures have already 
been formed; for this reason, their distribution will be 
correlated with the ISW signal, while the primary CMB 
anisotropy signal can't. It's for this reason that the best 
way to measure the ISW signal in the CMB anisotropics 
spectrum is the cross-correlation technique between the 
whole signal and a survey of the matter distribution in 
the universe |9| . The ISW temperature fluctuation in the 
direction n is given by 



f , s d$> 
Q ISW (n) = -2 / e~ T ^ — (n,z)dz, 



(12) 



where $ is the newtonian gauge gravitational potential 
and e -r ( z ' is the visibility function to account for a pos- 
sible suppression due to early reionization. Conversely, 
the observed overdensity in a given direction is 



(m) 



ip(z)5 m (rh, z)dz, 



(13) 



where S m is the matter density perturbation, b g the galac- 
tic bias and <p(z) the specific survey selection function. 
With these two definitions, one can compute the 2-points 
angular cross-correlation function 



JSW-LSS 



W = (Qisw(n)S LSS (m)), 



(14) 
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where ■d is the angle between the directions n and m, and 
the cross-correlation power spectrum |16| 



C, 



ISW-LSS 



/ f A> R l[ sw (k)l[< ss (k), (15) 



where An is the primordial power spectrum and 



tISW i 



{ft) = -2 I e- T(2) ^j ; [fcr(z)]dz (16) 



dz 



lt SS {k) = b g / vW^zMkrWdz, (17) 



where and are the Fourier components of the gravi- 
tational potential and matter perturbation, r(z) is the co- 
moving distance at redshift z and ji[kr(z)] are the spher- 
ical Bessel functions. 

In Figure 1 we plot some theoretical ISW-LSS power 
spectra computed under the case of generalized and silent 
chaplygin gases respectively. We consider the Chaply- 
ging gas as a candidate for dark energy alone, we fix 
Q m = 0.3, w = —0.9 and we let vary a. The prediction 
is made for an all-sky CMB survey plus a LSS survey 
centered at redshift = 0.3 with dispersion Az ^0.1. 
In the top panel we consider predictions for the Gener- 
alized Chaplygin. As we can see there is a strong, non 
trivial, dependence on a. Increasing a shiftss the epoch 
of matter to dark energy transition, affecting the ampli- 
tude of the ISW signal. At the same time, the Chaplygin 
gas affects the growth of perturbations in the CDM com- 
ponent, decreasing the growth in the case of a = 0. In 
the ISW — LSS cross signal the two mechanism are in 
competition and the signal is strongly dependent from 
the value of a and from the details ofthe LSS survey. 
On the bottom panel we plot the same predictions but 
in the case of the silent Chaplygin. As we can see, the 
additional intrinsic entropy perturbation has the effect of 
making the cross signal substantially independent from 
o. 
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FIG. 1: Predictions for the ISW-LSS power spectrum in the 
case of the generalized (top panel) and silent (bottom panel) 
Chaplygin. We fix £l m — 0.3 and w = —0.9. The LSS survey 
is centered at z» = 0.3 ± 0.1. 



III. CONSTRAINTS FROM DATA 

We perform a likelihood analysis using the collection 
of data presented in Gaztanaga et al. , which has the 
advantage of being publicly available and easy to imple- 
ment. In order to avoid degeneracies between the pa- 
rameters we consider only two possible viable Chaplygin 
dark energy models: a generalized Chaplygin model with 
a = and the silent quintessence model. For those mod- 
els we vary two parameters: the energy density Clchap 
and the equation of state w. 

The data we consider consist of measurements of the 
average angular cross-correlation around 9 = 5° between 
WMAP temperature anisotropy maps and several LSS 
surveys. The angular average around 9 = 5° ensures 
that the signal is not contaminated by foregrounds such 
as the SZ or lensing effects which are relevant at smaller 
angles (9 < 1°). Possible systematic contaminants, such 



as extinction effects, seem not to affect these data and 
for a more detailed discussion we refer to |7j. The data 
span a range of redshift 0.1 < z < 1 and for each redshift 
bin the data include an estimate of the galaxy bias with 
20% errors. These biases are inferred by comparing the 
galaxy-galaxy correlation function of each experiments 
with the expectation of best fit model to WMAP power 
spectra. There is little sensitivity to the scalar spectral 
index n s , while the dependence on the baryon density Qb 
can be non-negligible. In fact the presence of baryons in- 
hibites the growth of CDM fluctuations between matter- 
radiation equality and photon-baryon decoupling causing 
the matter power spectrum to be suppressed on scales 
k > k eq for increasing values of f2t> (k eq is the scale which 
enters the horizon at the equality). Over the range of 
scales which contribute to the ISW-correlation (k ~ 0.01) 
the sensitivity on f2b is still present. In order to limit 
the number of likelihood parameters we therefore assume 
a Gaussian prior on the baryon density fib = 0.04 and 
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we take the Hubble parameter as h — 0.7. We also as- 
sume a scale invariant primordial spectrum n a = 1 and 
fix the optical depth to reionization to WMAP best fit 
value t = 0.17 (again the ISW is not particularly affected 
by a change in those parameters). We marginalize over 
the normalization amplitude A s , although we found no 
difference assuming the WMAP value. In fact chang- 
ing A s shifts the overall amplitude of the angular cross- 
correlation of the same amount over different redshifts 
but it does not change the redshift dependence of the 
signal. Since the experimental data are corrected for the 
bias by comparing the measured galaxy-galaxy correla- 
tion function in each redshift bin to the WMAP best fit 
model, we rescale these biases to each of the dark energy 
model in our database as described in Q and Q . 

We compute for each theoretical model the angular 
cross-correlation as described in the previous section us- 
ing the selection function 



z 2 exp | 



-(z/z, 



\1.5l 



(18) 



where Zi is the median redshift of the i-th survey. Then 
following , we compute the average cross-correlation in 
the i-th bin, Cf , around 9 ~ 5°. 

The data points are an average over angles and are in- 
fered from surveys whose selection functions may overlap 
in redshift space, hence they are not independent mea- 
surements and indeed are affected by a certain degree of 
correlation. Since we have no access to the raw data we 
have no way of accounting for the first type of correlation, 
while using Eq. (|18|l we can estimate the correlation be- 
tween different redshift bins. We compute the correlation 
matrix p = {pij}, where pij is the fraction of overlapping 
volume between the i-th and j-th surveys (i.e. the diago- 
nal components are pu = 1). To be more conservative we 
have assumed that the different surveys cover the same 
fraction of sky, in general this is not the case and the frac- 
tion of overlaping volume can be smaller. We found that 
only two data points are highly correlated, since their se- 
lection functions overlap for about 70% in redshift space, 
while the correlation among the remaining data points 
are less than 20%. 

We compute a likelihood function C defined as 
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FIG. 2: Likelihood contour plots in the Q c -w plane for a 
generalized Chaplygin gas with a — (Top Panel) and for 
the silent Chaplygin case (Bottom Panel). The contours are 
68% (dotted), 95% (short-dashed), 99% (long-dashed) c.L. 



where Cf are the data and M^ 1 = pij/(ai<Tj) is our 
estimate of the inverse of the covariance matrix, which 
takes in to account the experimental errors and possi- 
ble correlation amongst the data, with <Ji the measured 
uncertainty in the i-th bin. 

The results of our analysis are plotted in Fig. 2 As 
we can see, in the case of the generalized Chaplygin with 
a = the current data constrain f2 c > 0.55 and w < —0.6 
at 95% c.L. Considering the silent Chaplygin model we 
obtain similar results with f2 c > 0.55 and w < —0.65 
again at 95% c.L. 



IV. CONCLUSIONS 

In this paper we have investigated the possibility of 
constraining Chaplygin models with the current Inte- 
grated Sachs Wolfe effect data. In the case of a fiat 
universe we found that cosmological viable generalized 
Chaplygin gas models must have an energy density such 
that f2 c > 0.55 and an equation of state w < —0.6 
at 95% c.L. We also investigate the recently proposed 
Silent Chaplygin models, constraining f2 c > 0.55 and 
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w < —0.65 at 95% c.L. These constraints have been 
obtained with a prior on the value of the baryon density 
ilb in agreement with BBN nucleosynthesis. Better mea- 
surements of the CMB-LSS correlation will be possible 
with the next generation of deep redshift surveys. For 
an ideal whole sky experiment the ISW signal may be 
detected at the level of ~ 10o\ This would translate in a 
factor 2/3 better determination of the results presented 
here and will provide independent and complementary 



constraints on unified dark energy models such as the 
Chaplygin gas. 
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